Phase diagram for Coulomb-frustrated phase separation in systems with negative 

short-range compressibility 
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Using numerical techniques and asymptotic expansions we obtain the phase diagram of a paradig- 
matic model of Coulomb frustrated phase separation in systems with negative short-range compress- 
ibility. The transition from the homogeneous phase to the inhomogeneous phase is generically first 
order in isotropic three-dimensional systems except for a critical point. Close to the critical point, 
inhomogeneities are predicted to form a BCC lattice with subsequent transitions to a triangular 
lattice of rods and a layered structure. Inclusion of a strong anisotropy allows for second- and 
first-order transition lines joined by a tricritical point. 

PACS numbers: 71.10.Hf, 64.75.-g, 64.75. Jk 



The appearance of spatial inhomogeneities is a gen- 
eral phenomenon occurring in a wide variety of sys- 
tems with competing interactions on different length 
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171 |l8(. In electronic systems strong electron correla- 
tions often drive the electronic compressibility negative 
as in the Hubbard model in infinite dimensi on [3J , the Fal- 
icov Kimball model[10(, the electron gas lift |20| , cuprates 
models 0, IE S], manganite mo dels [7L IE Inf. semiconduc- 
tor heterostructures|8(, to mention a few|2|. This signals 
a tendency to phase separation at macroscopic lengths 
which is frustrated by the long-range Coulomb inter- 
action. Although there is agreement that the system 
can become inhomogeneous at a mesoscopic scale if the 
Coulomb frustration is not too large, the nature of the 
transition has not yet been settled. This is an important 
question in general but particularly so in metallic systems 
where the existence of a second order quantum critical 
point separating an homogeneous phase from an inho- 
mogeneous one can disrupt the Fermi liquid behavior 

In this work we study a 4 model augmented with long- 
range interactions as a generic model of Coulomb frus- 
trated phase separation in systems with negative short- 
range compressibility. The model (or closely related vari- 
ants) has been used to describe inhomogeneities in a va- 
riety of systems [IE 

E1E1EEEE E3] including mixtures 

of block copolymers [131], charged colloids in polymeric 
solutions [l3| and electronic systems [13, EE] . We deter- 
mine the transition line from the homogeneous state to 
the inhomogeneous state using numerical and analyti- 
cal techniques and show that generically the transition is 
first order in three spatial dimensions (3D) except for a 
critical point (Fig. [T|) . This outcome changes if a strong 
anisotropy is taken into account. Then both first and 
second order transitions are allowed separated by a tri- 
critical point (Fig. [3]). In addition we study the crossover 
from harmonic to anharmonic inhomogeneities and the 
different topological transitions of the 3D isotropic sys- 



FIG. 1: (Color online) Phase diagram in three-dimensional 
isotropic systems. The small dots indicate the Gaussian 
instability line Q g . The thin (thicks) lines represent first- 
order transitions in the strong (weak) coupling approxima- 
tion. In the two limits they overlap with the corresponding 
numerically determined transition lines from homogeneous to 
droplet-like inhomogeneities (□), from droplets to rods (A), 
and from rods to layers (<()). 



tern. Close to the critical point the inhomogeneities are 
shown to form a BCC lattice. 

The model is defined by the following Hamiltonian 



H = J dx[0(x) 2 -l]' 
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with the scalar classical field 4> representing the local 
charge density, <f> the average density. A rigid background 
ensures charge neutrality. In the case of electronic sys- 
tems with a negative compressibility the model describes 
phenomena at large length scales compared to the under- 
lying lattice constant. It can be derived by expanding the 
coarse grained energy of the system around a reference 
density belonging to the negative compressibility density 
region. Fcrmionic effects like the Pauli exclusion princi- 
ple arc reflected in the coefficients. The reference density 



2 



0.5 


-0.5 
-1 



' 0.99 1 


\ 

- 


0.5 - - f , 




0.1 - / ' 




0.05 /' 


-(a) 


\ / 1 

V - / 









100 


Amplitude 


$5/ 


4>7 








50 












G/S £ 0.8 








-0.5 



0.5 



0.2 



FIG. 2: 



0.4 0.6 

Q/Q c 



0.8 





x/R c 

(Color online) (a) Behavior of the charge density 
modulation for <f> — and different couplings (labeled by the 
value olQ/Q c ). Near the Gaussian instability the modulation 
is close to a SCDW (thin full line). As Q decreases unhar- 
monicity is build in and the charge modulation tends to a 
square wave which has as limiting case Maxwell construction 
at Q = 0. (b) Evolution of the periodicity of the charge den- 
sity wave measured in units of £. Inset: Evolution of the 
harmonic amplitudes <f> n / X^ n '>o c t >n ' as a function of Q. The 
dots indicate the amplitudes for a layered structure of (j> — ±1 
as ruled by Maxwell construction at Q = 0. Amplitudes for 
even n vanish by symmetry. 



can be fine tuned so as to eliminate a small cubic term. 
(Our results can be trivially extended to take into ac- 
count a large cubic term) . Constant and linear terms are 
inessential, thus one obtains a double well energy which 
can be taken to the dimensionless form of Eq. ([T]) by mea- 
suring energies in terms of the barrier height, distances 
in terms of the bare correlation length £, and density <j> 
in units such that the double well minima are at </) = ±1. 
This leads to a renormalized Coulomb coupling in 3D, 
Q 2 = 4e 2 £ 2 /(eo|a|). Here e is the electron charge, e is 
a dielectric constant due to external degrees of freedom 
and a < is proportional to the inverse short-range com- 
pressibility. More precisely, it is the second derivative of 
the short range part of the energy per unit volume with 
respect to the density at the reference density. 

In general a and £ will depend on external parameters 
like pressure. They can even be taken as temperature 
dependent, as in Landau theory, in which case Ti. has to 
be interpreted as a free energy and the model becomes a 
mean-field description of a temperature driven transition 
to an inhomogeneous state. This can be useful, for ex- 
ample, to model inhomogcncities appearing below some 
temperature in manganites[2lj]. 

By computing the static response to an external field 
in momentum space, we get the charge susceptibility at 
finite momentum k which depends only on the modulus 
of the momentum, k, measured in units of 



X(k) 
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The susceptibility has a maximum at ko = [27r<3 2 ] 
which diverges as Q approaches the Gaussian instabil- 
ity line Q g (the dotted line in Fig. [T]) from above, where 
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and Q c = 1/ V27T- This indicates an 
instability of the homogeneous phase toward a sinusoidal 



charge density wave (SCDW) of periodicity 2R C = 2n/ko 
with vanishing wave amplitude at the transition and di- 
rection chosen by spontaneous symmetry breaking. 

The Gaussian transition can not survive at low Q. In- 
deed as Q — ► the Gaussian theory predicts inhomo- 
geneities in the range \<f>\ < l/Vo as opposed to Maxwell 
construction at Q = which predicts a globally inhomo- 
geneous state in the range \<j)\ < 1. We now show that 
the system never reaches the Gaussian instability except 
for a critical point in the phase diagram: the Gaussian 
line is preempted by a first order transition. Restricting 
to periodic textures, the free energy difference respect to 
the uniform state can be written in Fourier space as: 
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where the G's are the reciprocal lattice vectors and V is 
the unit cell volume of the periodic solution. 

The appearance of a self-generated cubic term in 
Eq. ([3]) calls for the possibility of a first-order tran- 
sition which can be treated analogously to the liquid- 
solid transition 22 , 23]. Assuming that the instability 
is weakly first order, which can be checked a posteriori, 
we restrict the sum to wave- vectors with magnitude G 
close to ko- To have an energetic advantage from the 
cubic term of Eq. ([3]) we need to find reciprocal lattices 
with triads of wave-vectors forming equilateral triangles 
so that they add to zero. By requiring symmetric struc- 
tures with inversion only three sets of vectors Lq are 
allowed that correspond to FCC, planar hexagonal, and 
icosahedral reciprocal lattices [23| . The free energy reads: 



— = (G)mcj) 2 G + 8(j)pm(j) 3 G + (j)c [3m [m- 1) + i 



6qm] 
(4) 

where the Fourier component amplitudes <\> G depends 
only on the modulus of G, m is the number of vectors 
in Lq and p (q) is the number of triangles (non-planar 
diamonds) to which each vector belongs when the tri- 
angles are accommodated in regular geometrical objects 
[24(. Upon minimizing Eq. (j4]) with respect to the wave- 
amplitude 4>g and G one finds G = kg. Equating the 
energies for the different structures one finds three first 
order transition lines between phases X and Y of the 
form 
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where X,Y = U, B, T, L stand for uniform, BCC, tri- 
angular and layered structures respectively (see below). 
The three first order transition lines [thin lines in Fig. [T] 
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join at the critical point (<j>, Q) = (0, Q c ) shown with a •. 
At the critical point one recovers a SCDW second-order 
phase transition with a charge susceptibility divergence 
since the cubic term of Eq. ([3]) vanishes. Away from the 
critical point the Gaussian line is the limit of metasta- 
bility of the homogeneous phase (doted line) . For all the 
phases, which we describe next, the order parameter at 
the transition goes linearly to zero as 4> —> confirm- 
ing that the transition is weakly first order close to the 
critical point. 

Approaching the first order lines from above the first 
inhomogeneous structure to become stable corresponds 
to the FCC reciprocal lattice defined by the m = 12 wave- 
vectors GV2/G = (±1,±1,0), (±1,0, ±1), (0,±1,±1) 
for which p = q = 2. This corresponds to a BCC crystal 
of inhomogeneities in real space with a.u,B — 103/45 rs 
2.29. 

Decreasing Q at finite <p the planar hexagonal lattice, 
with m = 6, p = 1, q = and asr ~ 3.44, becomes 
stable corresponding to rod-like inhomogeneities forming 
a triangular lattice similar to an Abrikosov lattice in a 
type II superconductor (middle thin line in Fig. [I]). For 
weaker Coulomb couplings and close to <j) = we find 
a subsequent morphological transition (lower thin line) 
that restores the translational symmetry in an additional 
direction and leads to a layered structure (m = 2, p — 
q = 0) with a T , L = 87/(19 - 6-y/6) 20.22. We also find 
that the icosahedral reciprocal lattice which corresponds 
to an icosahedral quasicrystal never becomes favorable. 

The above transitions [c.f. Eq. ((SJ)] are expected to be 
asymptotically exact close to the critical point. Instead 
for Q — > they are clearly inaccurate since, still, the 
range of stability of the inhomogeneous state is smaller 
than Maxwell construction. In this limit the weak first- 
order character of the transition is lost and more har- 
monics should be taken into account. 

A good approximation for small Q consists in assum- 
ing domains of uniform density of one or the other phase 

separated by sharp interfaces [2, 12, ee mm. tws IS 

a reasonable approximation for Q — > since there is a 
strong separation (Is >> Id >> £) between the typical 
interface scale length £, the typical size of the domains 
l d = (,/Q 2 / 3 and the screening length Is = Z/i^^Q) 
which controls the relaxation of the charge inside a do- 
main. The transition lines in this approximation are 
shown with the thick lines at the bottom of Fig. [TJ Again 
one finds the same topological transitions as in strong 
coupling but now the inhomogeneities form sharply de- 
fined spherical drops, cylinders and layers. The lattice is 
treated in the Wigner-Seitz approximation (WSA)[EEil- 
In order to study the crossover from the weak to strong 
coupling we have numerically minimized a discretized 
version of the energy, also in the WSA. For rod-like 
and droplet-like inhomogeneities we assume respectively 
cylindrical and spherical symmetry in order to reduce 
the minimization procedure to a one-dimensional effec- 
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FIG. 3: (Color online) The phase diagram for anisotropic 3D 
systems. The thin line corresponding to the Gaussian insta- 
bility determines the second order transition line above the 
tricritical point (•). Below the tricritical point the transition 
is first order and is determined numerically (□) and in the 
weak coupling limit assuming sharp interfaces @, H3] (thick 
line). Below the tricritical point the Gaussian line becomes 
the limit of metastability of the uniform phase. 



tive problem. The first-order transition points numeri- 
cally located are shown with the □, A, in Fig. [1] Wc 
find that as the coupling is decreased the size of the do- 
mains become much smaller than their distance which 
should make the approximation particularly accurate at 
weak coupling. Despite of that, the numerical result con- 
verges to the asymptotic expansions both at weak and 
strong coupling indicating a range of applicability of the 
WSA wider than expected. 

If one relaxes the WSA, the BCC lattice, found at 
strong coupling, is expected to evolve into a BCC Wigner 
crystal of drops at weak cou plin g which is the lowest 
energy Wigner crystal lattice Tj^ . Other lattices how- 
ever (FCC, HCP) are very close in energy which suggests 
that amorphous configurations will be very competitive 
as well. 

It is interesting to see how unharmonicity is build into 
the system as the coupling is decreased. In Fig. [2] (a) we 
show the charge profile for <j> = and different couplings. 
The SCDW smoothly evolves into the domain morphol- 
ogy that has as a limiting case the macroscopically phase 
separated state at Q — 0. Notice that the horizontal axis 
is normalized by the cell periodicity 2i? c , which becomes 
of the order of the linear size of the system as Q — ► 
[Fig- 0(b)]. 

The inset of Fig. 0(b) shows the behavior of the Fourier 
components. Close to the second-order phase transition 
the order n harmonics, 4> n = 4>G n with G n = nir/R c , 
behave as <f> n oc 

(Qg - Q) |n/21 . This follows from the 
fact that the modulation <p n couples with (</>i)'™' in the 
quartic term of Eq. (|3|). Higher harmonics proliferate 
as Q is decreased and converge to a rectangular profile 
corresponding to macroscopic phase separation at Q = 
(shown by • in the inset). 

The phase diagram changes dramatically if the gradi- 
ent term is made anisotropic i.e. £||/£_L > 1 where £m 
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(£_l) is the bare correlation length in the "hard" ("soft") 
direction. This can originate from an underlying crys- 
tal which favors certain orientations of the interfaces. In 
Fig. [3]we show the 3D phase diagram with two hard direc- 
tions and one soft direction for — * °o so that only 
one-dimensional modulations are allowed. The units are 
fixed as before with £j_ replacing £. 

In this case the cubic term in Eq. |3]) has no effect and 
one recovers the Gaussian instability where SCDW ap- 
pears at a second order transition in an extended range 
of coupling. Thus, for a metallic system, Landau damp- 
ing becomes relevant and one should take the coupling 
to the fermions explicitly into account 

For the same argument as before the Gaussian line 
can not persist up to Q — 0. We find, in fact a tri- 
critical point at Gf>, Q) = (±\/3/5, 16/25Q C ) where the 
transition becomes first-order (• in Fig. [3]). The position 
of the tricritical points and the behavior of the charge 
density modulation around it, can be studied retaining 
only two-harmonics in the order parameter with collinear 
wave-vectors : 

(x) — 4> + 2[4>i cos (Gix) + 4>2 cos (G 2 x)] 

and assuming (f>2 << 4>\- 

By expanding the energy around the tricritical point 
and minimizing with respect to the second-harmonic am- 
plitude <p2 one obtains an effective Landau free energy 
expansion for the modulated phase in terms of </>i only, 

£^I = r^+«40i + U6 ^ (6) 

with r = 4(Q~Q g ) /Q c , u 4 = 6-320 2 Q c /Q > and u 6 
a positive constant. This free energy has the canonical 
form for a tricritical point which as usual is determined 
by the vanishing of the quartic term on the Gaussian 
line@. 

In the first-order transition region [u^ < 0), the ap- 
pearance of the inhomogeneities is determined by the co- 
efficient uq in the free energy expansion, Eq. ©, whose 
precise value depends upon higher-order harmonic con- 
tributions. As before we have determined the first-order 
transition line by numerical minimization (□ in Fig. [3]) 
and from the weak coupling expansion assuming sharp 
interfaces @, H3| (dashed line). 

In conclusion, we have obtained the phase diagram 
and characterized the crossover of inhomogeneities from 
strong to weak coupling in a model of Coulomb frus- 
trated phase separation relevant for systems with a nega- 
tive short-range compressibility. For isotropic systems we 
find that the transition from the uniform phase to the in- 
homogeneous phase is always first order except for a crit- 
ical point. Close to the critical point inhomogeneities are 
predicted to form a BCC lattice with a subsequent tran- 
sition to a triangular lattice of rods and finally to a lay- 
ered structure. The transition lines continuously evolve 



into the weak coupling limit. In the case of anisotropic 
systems the transition to the inhomogeneous state can 
become second order with a tricritical point separating 
the second order line at strong coupling from the first or- 
der line at weak coupling. Inclusion of an explicit cubic 
term in Eq. (fT]), not considered here, will make the phase 
diagram asymmetric maintaining the topology. Thus we 
expect our results to be qualitative valid for a wide range 
of systems with Coulomb frustrated phase separation. 

One of the authors (CO.) acknowledges M. Beccaria 
for useful discussions. 
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